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The experimental observation that vanadate superlattices (LaVC>3) m /SrV03 show ferromag- 
netism up to room temperature [U. Liiders et ai, Phys. Rev. B 80, 241102R (2009)] is investi- 
gated by means of density functional theory. First, the influence of the density functional on the 
electronic and magnetic structure of bulk LaV03 is discussed. Second, the band structure of a 
(LaV03) m /SrVC>3 slab for m = 5 and 6 is calculated. Very different behaviors for odd and even 
values of m are found: In the odd case lattice relaxation results into a buckling of the interface 
VO2 layers that leads to spin-polarized interfaces. In the even case a decoupling of the interface 
VO2 layers from the LaO layers is obtained, confining the interface electrons into a two-dimensional 
electron gas. The orbital reconstruction at the interface due to the lattice relaxation is discussed. 
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I. INTRODUCTION 

The combination of low dimensionality and strong 
correlations proved fruitful in the quest for new ma- 
terials with functional-oriented properties. It lead to 
the high-T c superconductivity*^ transparent conduct- 
ing oxides* 3 - high capacitance heterostructures^ and large 
thermopower^ to quote a few. Recently, the creation of 
high-mobility electron gases with reduced dimensionality 
has also been pursued in artificially layered oxide sys- 
tems. They include (i) the interface between the two in- 
sulators SrTi0 3 and LaA10 3 ,— ~— (ii) systems with ultra- 
thin doping layers as Nb : SrTi0 3 in SrTi0 3 £ and RO 
layers with R = La, Pr, Nd in SrTi0 3 ^ and (hi) the 
surface of SrTi0 3 pii In some of these systems functional 
properties were found, such as superconductivit}^ and 
ferromagnetismJ^ 

Charge carriers with reduced dimensionality can be 
introduced in other complex oxides, too. For instance, 
ultra-thin SrV0 3 layers can be intercalated in LaV0 3 ^ 
resulting in ferromagnetism up to room temperature^ 
although none of the two oxides is magnetic at this tem- 
perature. This points out the relevance of the layered ge- 
ometry, since the solid solution Lai_ x Sr T V0 3 shows the 
typical phase diagram of a doped Mott-insulator with 
a transition from an antiferromagnetic insulator to a 
non-magnetic metal at about x = 0.2i 16 i 17 Contrary to 
the solid solution, (LaV0 3 ) m /SrV0 3 superlattices show 
ferromagnetic behavior above room temperature for m 
ranging from 3 to 6, with a higher magnetization for 
even m values than for odd values. Thus, the geomet- 
rically confined doping in the superlattices gives rise to 
a magnetic phase which is not observed in the randomly 
doped solid solution. 

From the theoretical point of view, the reduced di- 
mension of electronic states confined in a potential cre- 
ated by the interface of two insulators, or by geometri- 
cally confined doping, leads to a renewed interest in two- 
dimensional systems. Studies of the three band Hubbard 



model predict a ferromagnetic ground state in case of 2/3 
and 3/4 band filling for a reduced bandwidth due to the 
enhanced correlations of the electrons^ Ferromagnetism 
was also theoretically predicted for a model specific to 
two t2 9 orbitalsf^ at SrTi0 3 /LaV0 3 interfaces 2 ^ and in 
SrTi0 3 /SrV0 3 superlattices* 2 ^ indicating a strong cou- 
pling of the electronic, magnetic and orbital degrees of 
freedom to the geometry of the supcrlattice including the 
number of SrV0 3 layers. Again, in both cases the mag- 
netic phase was largely explained by the confinement and 
a change in the orbital occupations. 

However, in the case of superlattices or heterostruc- 
tures, the reduced bandwidth and therefore enhanced 
correlations are just one part of the story. As in those 
systems two different materials are epitaxially grown on 
top of each other, also the misfit of the lattice parame- 
ters of both materials plays a role for the properties of the 
system mediated by the distortion of the structure due 
to strain. Especially in oxides the distortion can have 
a strong influence on the properties since, as mentioned 
above, the electronic, spin and lattice degrees of freedom 
are strongly coupled. This effect was shown theoretically 
in distorted LaV0 3 where, depending on the distortion of 
the lattice, different antiferromagnetic structures become 
energetically favorable^ 2 - 

In order to analyze the ferromagnetic state observed in 
(LaV0 3 ) m /SrV0 3 superlattices, electronic band struc- 
ture calculations are helpful to confirm the confinement 
of the electronic states and to elucidate the origin of the 
ferromagnetism. For this purpose, we will use density 
functional theory and address the distortion of the lat- 
tice. Furthermore, the superlattices show an asymmetry 
of the magnetic properties between odd and even num- 
bers of LaV0 3 unit cells. Thus, one aim of the calcula- 
tions is to provide evidences for a possible difference in 
confinement between the two geometries. 

The paper is organized as follows. In Sec. [H] we discuss 
the influence of the density functional on the electronic 
and magnetic structure of bulk LaV0 3 . The structure 
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FIG. 1: (Color online) Top: Proposed spin and orbital order in the rare earth vanadates. Bottom: Density of states projected 
on the V 3d orbitals. Left: C-type spin ordered LaVOs. Right: G-type spin ordered LaV03. The orbital occupancy of the VI 
atom in the left lower corner (red line) is compared to the orbital occupancy of the according V of the upper layer with the 
same spin orientation (blue line). The latter DOS is multiplied by —1 and plotted along the negative y-axis. The c/a ratio is 
1.02. 



model of (LaV03) m /SrV03 superlattices is addressed in 
Sec. IIII1 followed by the electronic structure, the charge 
carrier distribution and the orbital occupancy of the m = 
5 and m — 6 superlattices, as representative of even and 
odd superlattices, respectively. Finally, our findings are 
summarized in Sec. IIVI 



II. ELECTRONIC AND MAGNETIC 
STRUCTURE OF BULK LaVOs 

There are many systems for which there is a consensus 
about the appropriate density functional. This is not the 
case for the vanadates and we therefore start our consid- 
erations addressing this point. The rare earth vanadates 
with perovskitc structure, including LaV03, have the V 
ion in a 3+ valence state, i.e., the two electrons in the 
V 3d shell occupy two of the tig orbitals with S = 1. 



The crystal field splitting of the 3d orbitals in threefold 
degenerate t2 g and twofold degenerate e g orbitals is a 
consequence of the octahedral coordination of the V ions 
by O atoms. The rare earth vanadates show simultane- 
ous antiferro-type spin and orbital ordering^ The spin 
ordering of LaV03 is C-type, i.e., the alignment is an- 
tiferromagnetic within the xy-plane and ferromagnetic 
along the z-direction^ The accompanying orbital order 
is G-type, i.e., it alternates in all directions. Other vana- 
dates show G-type spin ordering together with C-type 
orbital ordering. For an overview of the phase diagram 
see Ref. [24|. It has been shown in Ref. [22| that strain due 
to the coherent growth on a substrate can strongly af- 
fect the orbital ordering. For an increasing c/a ratio the 
G-type order is destabilized. 

It is essential that the approximation to the density 
functional treats the lattice, magnetic, and electronic 
properties on an equal level. While the local density 
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approximation (LDA) typically underestimates the bond 
lengths and magnetic couplings, the generalized gradi- 
ent approximation (GGA) overestimates both. For bulk 
LaVC>3, the GGA correctly predicts the C-type spin 
ordering, while the LDA tends to delocalize the elec- 
tronic states, which is critical in the case of narrow 3d 
bands. To obtain an insulating state, corrections be- 
yond the LDA/GGA have to be considered, such as the 
LDA/GGA+C/ approach in which a local interaction is 
added for the 3c? orbitals^ This enhances the separa- 
tion of the 3d energy levels and opens a band gap. It is 
found that the magnetic ordering can be changed by this 
treatment^ 

A further step towards a rigorous treatment can be 
achieved by hybrid functionals which incorporate a non- 
local Fock exchange. Hybrid functionals such as PBEO 
and HSE yield an electronic structure comparable to 
self-consistent GW calculations for much lower compu- 
tational cost. Even the local EECE approach has shown 
clear advantages over the LDA/GGA+L/ in describing 
the electronic structure in nickelate stripe phases^ How- 
ever, it is known that the Fock term (without correction 
of the correlation term) results in a substantial overesti- 
mation of the magnetic coupling and, therefore, in an un- 
reliable determination of the magnetic ordering, includ- 
ing phase transitions. A GW treatment leads to accurate 
electronic structures and magnetic phase transitions but 
is computationally too costly for superlattice systems. 
The best compromise therefore is the LDA/GGA+C/ ap- 
proach, which we apply in the following. Specifically, we 
employ the WIEN2k augmented plane wave plus local 
orbitals code^ We find that the energies and electronic 
states obtained for the fully localized and around mean 
field double counting corrections show the same behavior. 

The interaction parameter U is adjusted as follows. 
For U > 5 cV the C-type spin ordering is energetically 
favored over the G-type ordering. A band gap of 1.5 eV 
is obtained for U = 6.12 eV in the case of C-type spin or- 
dering, while the same gap is found in the case of G-type 
spin ordering for U = 4.8 eV, see Fig. 1. In addition, 
the C-type orbital ordering under G-type spin ordering 
depends only weakly on U and on the lattice strain. On 
the other hand, the G-type orbital ordering under C- 
type spin ordering is very sensitive to strain. Disorder 
or ferro-type orbital ordering is found for c/a > 1 (which 
is fulfilled in our superlattice system). Ferromagnetism 
and A-type antiferromagnetism result in higher energies 
than the C- and G-type spin orderings for all approxi- 
mations under study. Electronic densities of states for 
C-type and G-type spin ordered LaV03 for c/a = 1.02 
are shown in Fig. 1. While the G-type density of states 
(DOS) indicates very narrow bands, the bands are much 
broader for C-type spin ordering despite the higher U 
value. For C-type spin ordering, the d xz orbital, forming 
a double peak, is occupied, while the d xy / yz orbitals are 
partially occupied and therefore subject to orbital order- 
ing, in agreement with Ref. I22T In contrast, for G-type 
spin ordering, the d xy orbital is occupied, with a much 



less pronounced DOS structure, and the orbital ordering 
affects the d xz j yz orbitals. 

Within GGA the ground state of bulk SrV03 is antifer- 
romagnctic with \iy = 0.2 fig, but the energy difference 
to the non-magnetic solution is small, amounting to 0.2 
meV per V atom^ For the GGA+C/ method, the C-type 
antifcrromagnet is favored by 30 meV per V atom over 
the A-type antifcrromagnet and by 140 meV per V atom 
over the ferromagnet. A local hybrid approach with 25% 
Fock mixing would favor G-type spin ordering, since a 
C-type starting configuration changes to G-type during 
the self-consistency calculation. The onsite interaction 
does not open a band gap in bulk SrV03 but the com- 
pound stays metallic. Around the Fermi energy, we find 
the t2 g orbitals equally partially occupied, while the e g 
orbitals are empty. Even though the spin-orbit coupling 
has little impact on the ground state of LaV03 it must be 
included in the case of SrVOs^i Indeed it splits the t^g 
states and gives rise to a band gap of 0.61 eV. In this case 
C-type spin ordering (magnetization in the x-direction) 
comes along with filled d xy and half-filled d XZl yz orbitals, 
yielding ferro-type orbital ordering. 



III. (LaV0 3 )m/(SrV0 3 ) SUPERLATTICES 

The structure setup for the calculations of the superlat- 
tice starts from the tetragonal phase of LaVOs, neglect- 
ing an orthorhombic lowering of the symmetry due to ro- 
tations of the oxygen octahedra. In addition, ideal atom- 
ically sharp interfaces are studied. We put one SrV03 
layer between m = 5 or m = 6 LaV03 layers, representa- 
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FIG. 2: (Color online) Supercell setup for the 
(LaV03) m /SrV03 heterostructure for m = 5 (left) and 
m — 6 (right). La (Sr) atoms are represented by blue (green) 
spheres and V 3+ (V 2+ ) ions are represented by orange (pink) 
spheres. The grey squares indicate the oxygen octahedra. 
The oxygen atoms themselves are omitted, except for the 
bottom interface layer for m = 5 where black spheres are 
shown to highlight the buckling of the VO2 layer. In the 
bottom of the m = 6 structure the (101) plane of the 
tetragonal (orthorhombic) cell is indicated in black (red). 
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FIG. 3: (Color online) Density of states projected on the (left) V 3d xy and (right) V 3dxz, y 
The inversion layer is denoted as bulk layer. The labeling of the layers is shown in Fig. [2] 
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tive of the odd and even superlattices, respectively. The 
lattice constants are a = 3.88 A and c = (m + 1)3.95 A. 

As observed experimentally, the c/a ratio in the super- 
lattice deviates from the bulk ratio of LaV03.— In bulk 
LaVC>3 it is smaller than 1, c/a = 0.98^ but the strain 
exerted by the SrVC>3 layers distorts the LaVC>3 sublay- 
ers in the superlattice. In case of the odd superlattice, 
six VO2 layers lie in between the SrO layers and the LaO 
forms the inversion layer, while in the even superlattice 
with its seven VO2 layers, the central VO2 layer is the 
inversion layer (sec Figure [2]). 

In a first step, we neglect the antifcrromagnetic order 
in the a6-planc and restrict ourselves to a 1 x 1 supercell 
in the afo-plane. In c-direction, ferromagnetic alignment 
is proposed, as expected for C-type spin order. It turns 
out in a second step that when taking into account also 
the magnetic ordering in the a6-plane by a 2 x 2 supercell 
the results concerning the interface electronic structure 
do not depend on the ordering in the a6-planc. 



A. Superlattice m = 6 

Let us start our investigations with the m = 6 super- 
lattice and first focus on the orbital occupancy recon- 
struction at the interface. Deep in the LaV03 region a 
V 3+ valency is realized, with the 3d xz and 3d yz orbitals 
occupied and the V 3d xy orbitals empty, as revealed by 
Fig. [3] On the contrary, for the interface layers (the two 
VO2 layers adjacent to the SrO layer) the V 3d xy orbital 
is occupied and the V 3d xz ^ yz orbitals are partially occu- 
pied, see Fig. [31 together with a small admixture of the e g 
orbitals. This results in a smaller valency of the V ions, 
which shows that they mostly retain the oxidation state 
of the respective bulk material. This provides us with a 
first indication that the doping of LaV03 is geometrically 



confined. 

This statement is further confirmed by the band struc- 
ture shown in Fig. |4j It exhibits four interface-based 
bands crossing the Fermi energy, two of them giving 
rise to fairly symmetrical Fermi surfaces (in the k x k y - 
plane) centered around the T point. These bands are 
characterized by rather large Fermi velocities. In addi- 

1-5 -1 1 . . ■ 1 . ; 1 




FIG. 4: Band structure of the m = 6 superlattice with non- 
magnetic interfaces. The interface V 2>d xz +3d yz orbitals are 
highlighted. 
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tion, small hole pockets centered around the M-point are 
found. They are characterized by clearly smaller Fermi 
velocities and are due to bulk states. Therefore, the fast 
charge carriers are mostly confined to the interface lay- 
ers, though in states also containing a small admixture 
from the sub-interface layer. Additional features of the 
band structure point to an interesting two-dimensional 
character of the bands crossing the Fermi energy. Prin- 
cipally constituted by d xz (d yz ) orbitals which disperse 
along the T-X (r~Y) direction, they unexpectedly show 
considerable dispersion along X-M (Y-M), pointing to 
an admixture of other orbitals to the d xz (d yz ) orbitals. 
Also, the dispersion around the T point is rather flat and 
not cosine-like as one would be expected for pure d xz /d yz 
bands. 

A partial occupation of the d xz and d yz orbitals es- 
tablishes the prerequisite of an orbitally ordered state. 
However, this is neither observed in our calculations nor 
in the experiments, where it was shown that these sys- 
tems resemble rather good metals at low temperature.— 
An explanation for this behavior may be found in the ob- 
served structure relaxations, where the most important 
ones are summarized in table HI In the interface layers the 
V (O) ions move towards (off) the SrO layers, resulting 
in a sizeable buckling of the interface layers with an am- 
plitude close to 0.11 A. The first LaO layers move quite 
homogeneously towards the SrO layer, with a buckling 
amplitude smaller than 0.01 A, while the buckling am- 
plitude of the sub-interface VO2 layer is close to 0.04 A. 

As a result, the bond lengths are strongly reduced be- 
tween the V ions in the interface VO2 layers and the O 
ions in the SrO layer, becoming smaller than a/2, i.e., 
the distance of the V ion to the O ions in the xy-plane 
of the sample positioned at a/2 is larger than that to 
the O ion in the z-direction. Consequently, the non- 
degenerate V energy level e xy lies below the e xz ^ yz levels 
in the interface layer, leading to an enhancement of the 
dispersion in those layers. On the other hand, the en- 
larged bond length between the interface layer and the 
adjacent LaO layer suppresses the orbital overlap in the 
z-direction and favors therefore an electronic decoupling 
of the interface VO2 layer with respect to the other VO2 
layers, as shown in Figure [SJ The same tendency of the 
buckling can still be seen in the sub-interface layer, but 



TABLE I: Displacements (in A) off SrO layers of the V and 
O atoms in the interface and middle/sub-interface layers for 
the heterostructures with m — 5 and m — 6. In the first LaO 
layers for m = 5 (m = 6) the motion towards the SrO layer 
of the La ions amounts to 0.118 A (0.093 A), and that of the 
O ions to 0.116 A (0.010 A). 
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FIG. 5: (Color online) V-0 and O-V bond lengths along the 
z-direction for m — 5 (top) and m = 6 (bottom). The vertical 
dashed line indicates a/2. 



with a much smaller amplitude. The O-V bond length in 
this case is near to c/2 as for the other VO2 layers. Thus, 
the larger distance between the V ions and the surround- 
ing O atoms is to be found along the ^-direction, inducing 
the bulk orbital occupancy and, hence, an orbitally or- 
dered insulating state like in the bulk. 

Thus, the calculated band structure accounts for a 
two-dimensional metallic state, but it cannot account 
for the experimentally observed ferromagnetism. In a 
Stoner picture this would be understood as the result 
of a rather small DOS at the Fermi energy. Therefore, 
should this system be ferromagnetic then ferromagnetism 
would likely be caused by further correlation effects, as 
taken into account in calculations of multi-orbital Hub- 
bard models. 18 ' 19 Moreover, we have verified that the 
above electronic structure is stable against the employed 
U value and lattice relaxation, as is the establishment of 
an ultrathin two-dimensional electron gas. 



B. Superlattice m = 5 

The m = 5 superlattice hosts different physics. Again, 
lattice relaxation plays an important role. The obtained 
buckling of the interface layers is summarized in Table HI 
The interface VO2 layer shows a buckling with the same 
displacement directions as for m — 6 (the V ions move 
towards the SrO layer and the O ions towards the LaO 
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layer) but the amplitude is smaller. Notably, as shown 
in Fig. [5j while for the m = 6 case the distance between 
the V ion in the interface layer and the O ion in the 
SrO layer becomes smaller than a/2, inducing the inver- 
sion of the orbital occupancy, in the m = 5 case the 
bond lengths stay larger than a/2 in all the VO2 lay- 
ers. The smaller amplitude of the buckling therefore does 
not induce an inversion of the orbital occupation and the 
confinement of the charge carriers is clearly less effective 
than in the even case. Furthermore, the reduction of the 
bond lengths reaches further into the LaVC>3 layer, where 
the bond length of the V ion with the second LaO layer 
counted from the SrO layer is still strongly reduced. In 
comparison, in the m — 6 case the same bond length is 
larger than c/2. 

The change of the structural relaxation is visible in 
the resulting layer-projected densities of states in Fig. [5] 
They consist of the spin-majority density of states, as 
no sizeable contribution from the spin-minority bands is 
obtained. All layers are clearly off stoichiometry and all 
V ions are in a mixed valency state, yet closer to the 3+ 
valency in the bulk layers than in the interface layers. 
Moreover, all V 3d electrons occupy the 3d xz and 3d yz 
orbitals, while all 3d xy orbitals are empty, in contrast to 
the m = 6 superlattice. 

The band structure of the m = 5 superlattice is shown 
in Fig.[7J As there are no spin- minority bands close to the 
Fermi level, only spin-majority bands are shown. There 
are two narrow bands formed by the interface V 3d xz ^ yz 
orbitals around the Fermi level. The four next lower ly- 
ing bands are due to the middle and bulk layers with 
a larger band width. The bands at the Fermi level are 
separated by a band gap of about 1 eV (measured at 
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the r point) from the parabolic bands of the bulk states. 
In addition, the bands with the largest Fermi velocities 
mostly trace back to the middle layer V 3d xz , yz orbitals, 
with some admixtures of the bulk and interface layers. 
Therefore, transport is no longer confined to the inter- 
face layers but there is a distinct modulation associated 
to it. There are fast charge carriers mostly confined to the 
next-to-intcrface layer, with a small admixture of their 
two neighbors, and slow charge carriers mostly confined 
to the interface layer. 

Narrow d bands with high DOS at the Fermi level of- 
ten explain ferromagnctism in the 3c? transition metal 
oxides. Here the high DOS follows from the slow charge 
carriers that are confined to the interface layers. We 
therefore have performed calculations for 2x2 supercells 
in the a6-plane and find that the ground state shows a 
ferromagnetic alignment of the spins within the interface 
layer, whereas they are aligned antifcrromagnctically in 
the other layers. 



IV. CONCLUSION AND OUTLOOK 

Our calculations show that the doping resulting 
from the charge carriers donated by the SrO layers 
in (LaV03) m /SrV03 supcrlattices leads to rich physics 
built on fast and slow charge carriers together with lat- 
tice relaxation. For both types of superlattices the lattice 
relaxation plays an important role, demonstrated by the 
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FIG. 6: (Color online) Spin-majority density of states pro- 
jected on the V 3d xz orbital. The labeling of the layers is 
shown in Fig. [2] 



FIG. 7: Band structure of the m — 5 superlattice with ferro- 
magnetically polarized interfaces. The interface V 3d xz +3d yz 
orbitals are highlighted. 
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obtained buckling of the VO2 interface layers and alter- 
nation of short and long V-0 bond lengths along the 
c-axis. Although the buckling of the interface VO2 lay- 
ers is qualitatively the same in the two studied cases, 
the critical difference in the amplitude leads to two dif- 
ferent band structures. In the odd m = 5 case, slow 
charge carriers are mostly confined to the interface layers 
in which they order fcrromagnctically. The confinement 
of the fast charge carriers is less effective and their distri- 
bution peaks in the middle layers. For even superlattices 
the relaxation confines the fast charge carriers to largely 
decoupled interface layers in which they form an ultra- 
thin two-dimensional electron gas. No magnetic order is 
found and the distribution of slow charge carriers peaks 
in the sub-interface layers. Therefore, coming back to 
our original question, it seems unlikely that the geomet- 
ric confinement of the charge carriers alone is the origin 
of the experimentally observed ferromagnetic state. In- 
stead, the fact that the even and odd superlattices show a 
different behavior points out the relevance of the number 
of unit cells of the LaVC>3 layers. 

In our calculations we did not take into account the or- 
thorhombic structure of LaMO^r^ generated by the oc- 
tahedral tilts which were shown to be present in pure 
LaVC>3 thin films as well.— The related unit cell is larger 
than the tetragonal cell used in the calculations, in par- 
ticular along the tetragonal c-axis, where the out of phase 
rotation of the octahedra leads to a doubling of the unit 
cell, see Fig. [2j Therefore, in an even superlattice it is 
possible to establish the bulk structure including the oc- 
tahedral rotations in the LaVC>3 layer. On the other 
hand, in the odd superlattice the thickness of the LaVC>3 
layer corresponds to a half-integer number of the or- 
thorhombic unit cell, so that the bulk structure (and 
therefore the bulk physics) of the LaVC>3 layer cannot 
be established. The compensation by the structure of 
the lack of the energetically favorable out of phase ro- 
tation may lead to the small but critical differences of 
the lattice relaxations and therefore the orbital occupa- 
tions between the even and odd case. The combination 
of the individual reactions of the structural, orbital, and 



electronic degrees of freedom leads in the end to a ferro- 
magnetic ground state in the odd superlattices. 

Our calculations confirm therefore the experimentally 
observed magnetization oscillations^ although in the ex- 
periment the even superlattices were found to show a 
higher magnetization than the odd ones. This difference 
can be explained by the fact that the thickness of the 
LaV03 in the experimental study is a nominal one, where 
the announced thickness does not take into account any 
intermixing of La and Sr during the growth. It was shown 
that the interface between the LaV03 and SrV03 layers 
is abrupt but exhibits steps of typically one pcrovskite 
unit celli22 The steps lead to the presence of Sr atoms in 
the first LaO layer next to the SrO layer. Taking into 
account the growth direction, it can be argued that the 
thickness of the pure LaV03 layer will thus be one unit 
cell less than the nominal one, resulting in an inversion 
of the nominally odd and nominally even superlattices. 

In conclusion, two scenarios for the even and odd su- 
perlattice have been found. In both cases, the d xz and d yz 
orbitals form bands crossing the Fermi level. In the even 
superlattice, the orbital occupancy changes from the in- 
terface to the bulk layers, while we find an homogeneous 
orbital occupation in the odd superlattice. The room 
temperature magnetism is related to the incomplete for- 
mation of the orthorhombic unit cell of LaV03 in odd 
superlattices. A next step could be a comparison of the- 
oretical and experimental macroscopic transport proper- 
ties of the confined fast charge carriers. 
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